Turbulent flows for relativistic conformal fluids in 2+1 dimensions 



<N 

O 
<N 

> 

o 

o 

00 



(N 

> 
<N 
O 

o 

(N 



X: 



Federico Carrasco^'^, Luis Lehner^'^, Robert C. Myers^, 
Oscar Reula^, Ajay Singh^'^ 
^FaMAF-UNC, IFEG-CONICET, Ciudad Universitaria, 5000, Cordoba, Argentina. 
^Perimeter Institute for Theoretical Physics, Waterloo, Ontario N2L 2Y5, Canada. 
^Department of Physics, University of Guelph, Cuelph, Ontario NIC 2W1, Canada. 
Department of Physics & Astronomy and Cuelph- Waterloo Physics Institute, 
University of Waterloo, Waterloo, Ontario N2L 3C1, Canada 

We demonstrate that relativistic conformal hydrodynamics in 2+1 dimensions displays a turbu- 
lent behaviour which cascades energy to longer wavelengths on both flat and spherical manifolds. 
Our motivation for this study is to understand the implications for gravitational solutions through 
the AdS/CFT correspondence. The observed behaviour implies gravitational perturbations of the 
corresponding black brane/black hole spacetimes (for sufficiently large scales/temperatures) will 
display a similar cascade towards longer wavelengths. 



I. INTRODUCTION 

The AdS/CFT correspondence [l|, Q provides a re- 
markable framework for studying certain strongly cou- 
pled gauge theories in d dimensions by mapping to weakly 
coupled gravitational systems in d + 1 dimensions. Of 
particular interest is the relation between gauge the- 
ory plasmas and black hole geometries Here the 
correspondence was used to show that perturbed black 
brane geometries at the classical level have a dual de- 
scription in terms of fluid dynamics equations govern- 
ing long- wavelength perturbations about an equilibrium 
state. This fluid/gravity correspondence builds on the 
highly successful and ongoing program to calculate the 
near-equilibrium transport coefficients for strongly cou- 
pled plasmas using holographic techniques, e.g., [6|, 0- 

This duality has been exploited in a large number of 
works which exploit known gravitational behaviour to in- 
fer properties of diverse systems described by strongly 
coupled field theories in the large N limit. On the hydro- 
dynamical front, work in [8] describes how the onset of 
naked singularities can be tied to finite-time blowups in 
hydrodynamics. Additionally, numerical simulations are 
increasingly being exploited to understand isotropization 
and thermalization of systems starting far from equilib- 
rium, e.g., jol-fTsj. 

In all these works, the approach has been to under- 
stand the behaviour of relevant systems in the gravita- 
tional side of the duality and infer from it properties of 
the gauge theory dynamics. In the present work, we fol- 
low the opposite route, namely to study particular phe- 
nomena that might arise in the field theory to under- 
stand possibly unexpected phenomena on the gravita- 
tional side. 

Our starting point is a simple and well-known ob- 
servation about the behaviour of turbulence of Newto- 
nian fluids in two spatial dimensions. Specifically in this 
case, turbulent behaviour induces an energy cascade from 
short to long wavelengths e.g., [I6|. This contrasts to 
the standard cascade from long to short scales charac- 
terizing turbulence in three and higher dimensions. (H^ 



The obvious question is therefore whether this effect 
appears in the relativistic conformal fluids relevant for 
the AdS/CFT correspondence, and if so, what is the 
dual interpretation in the gravitational theory. A pos- 
itive answer would seem to distinguish AdS gravity in 
four dimensions from higher dimensions in a unique way. 
An affirmative answer is also particularly intriguing as 
the two known instabilities in AdS spacetimes, super- 
radiance and the recently found 'mildly-turbulent' be- 
haviour in , both induce energy cascades from low to 
high frequencies for all d > 3. 

To answer this question, which has also been raised 
earlier [Bj, we examine the behaviour of specific fluid flows 
which are dual to perturbations of Schwarzschild and 
Kerr black holes in AdS 4,. We study the problem for 
spherical horizons in global coordinates, as well as pla- 
nar horizons described by Poincare coordinates. Specifi- 
cally, we study the relativistic Euler equations with the 
particular equation of state corresponding to a confor- 
mal fluid on a fixed 5R x S*^ manifold for the former case 
while 5RxT^ for the latter. In each background, we exam- 
ine perturbations of the stationary configurations dual to 
the corresponding stationary black holes. We then exam- 
ine the onset of turbulence and its cascade behaviour in 
this setup, and compare it with the expected behaviour 
for the case of incompressible non-relativistic flows with 
'standard' equations of state. We find that, as in the 
Newtonian case, turbulence leads to perturbations cas- 
cading from shorter to longer wavelengths for relativistic 
conformal hydrodynamics in 2+1 dimensions. 

This work is organized as follows: in the remainder of 
this section we briefly review some well known aspects 
about turbulence in two and three spatial dimensions. 
Section II discusses the initial configurations considered, 
as well as details of our numerical implementation. Sec- 
tion III presents results for the cases considered for con- 
formal hydrodynamics on the sphere (the results on the 
torus, which are qualitatively similar, are included in an 
appendix). We discuss the consequences of the obtained 
behaviour in section IV. (Visualizations of the fluid flows 
studied here can be found in plsf). 
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A. Turbulence 

Turbulence is a ubiquitous property of fluid flows ob- 
served in nature [19]. Qualitatively, we might describe 
turbulence as a flow regime characterized by chaotic or 
stochastic behaviour. Most of our theoretical under- 
standing of turbulence comes from the study of non- 
relativist ic incompressible fluids. Certainly, while a full 
understanding of turbulence is not yet available, some 
robust results do exist. Namely, for the inviscid case in 
two spatial dimensions, a global regularity theorem has 
been proved, together with theorems about uniqueness 
and existence of solutions [20] implying no singularity of 
the velocity field can develop in a finite time. 

For the three-dimensional case it is not yet known 
whether the same holds true and resolving this issue con- 
stitutes a major open problem (see e.g., [21|). It is known 
however, that qualitatively two and three dimensional 
turbulent fluid flows exhibit profound differences arising 
from the existence of a key conserved quantity which has 
a radical effect in the fluid's turbulent behaviour. This 
quantity, dubbed enstrophy^ — see discussion around 
(p!8|) — has been argued to imply a very different cas- 
cade picture in two dimensions, as compared to the three- 
dimensional one [22]. Small scales will support a ^direct 
enstrophy cascade' towards smaller wavelengths, with all 
the enstrophy dissipation taking place on the shortest 
scales. The energy flux towards small scales will then 
be damped and the energy will be, instead, transferred 
to larger scales in an 'inverse energy cascaded This be- 
haviour is observed in nature, controlled experiments and 
numerical simulations. 

In a phenomenological theory of two-dimensional tur- 
bulence [22I, [23I , the existence of two inertial ranges was 
pointed out: a direct enstrophy cascade at small 
wavelengths, and an inverse energy cascade with spec- 
trum at larger scales, were predicted. 

Numerical simulations have shown the emergence of 
strong coherent vortex structures that dominates the flow 
after some time [24]. These vortices emerge as anoma- 
lous fluctuations at small scales, and have a lifetime much 
longer than their characteristic eddy turnover time (25| . 
As the dynamics continues, two such vortices might col- 
lide and usually merge if they rotate in the same direc- 
tion, forming a vortex of larger scale. In the process the 
energy is nearly conserved, so it acts as the mechanism 
of inverse energy cascade. The statistical distribution 
of vortices over scales leads to an energy spectrum 
which is much steeper than the originally expected 
and vortices were recognized as the fluctuations responsi- 
ble of intermittency and possible anomalous dimensions 

0. 

The emergence and dynamics of vortices still pose a 
number of difficult questions, and the relevance of the 
initial data and external forcing is not yet fully under- 
stood. An exciting possibility, of course, would be that 
holographic studies might shed new light on this problem 
from a fundamentally different point of view. 



II. EQUATIONS, RATIONALE AND 
NUMERICAL IMPLEMENTATION 

Our goal is to study the dynamics dual to a perturbed 
black hole in 3 + 1 dimensions from the global point of 
view, so we consider the field theory in 2 + 1 dimensions 
on the sphere (or a torus, see Appendix A). To do so, 
and in order to deal with a well posed problem, we re- 
strict to the zeroth order expansion of the theory, i.e., 
to the equations determined by the conservation of the 
stress energy tensor of a perfect fluid. One reason for 
this choice is that the inclusion of viscous terms would 
yield an acausal system of equations pTMsoj (for a recent 
discussion in the holographic context see e.g., [3l|).[55| 
Hence properly incorporating these effects poses a serious 
complication for the corresponding simulations. 

At first sight, working with perfect fluid equations 
would seem to be a severe limitation for our study and 
the conclusions that can be drawn from this work. How- 
ever, recall that hydrodynamics treats the conservation 
of a stress-energy tensor in a gradient expansion, e.g., 
[32, 33]. For the conformal case, no intrinsic scales appear 
in defining the fluid and hence the temperature naturally 
controls the equation of state and all of the non- vanishing 
transport parameters. Hence to the stress tensor for a 
conformal fluid, in d spacetime dimensions, takes the fol- 
lowing form 



du^u^y) 



^a,.+0{T-')] (1) 



where aab describes the first-order viscous contribution 
and a, /3 are dimensionless coefficients which character- 
ize the fluid. From this expression, we see that the vis- 
cous terms can be arbitrarily suppressed by increasing 
the temperature. To be more precise, if the characteris- 
tic scales L controlling the flow are kept fixed, i.e., the 
size of the sphere (or torus) in the present case, the gra- 
dient expansion becomes an expansion in 1/LT and the 
higher order terms are suppressed by setting LT ^ 1. 
Thus the perfect fluid can be seen as an arbitrarily good 
approximation in this regime. Further, as we discussed 
above, turbulence is expected to generate a cascade from 
short to long scales for fluids in two spatial dimensions. 
In fact, our simulations will confirm this expectation and 
so the interesting dynamics here indeed progresses mainly 
towards longer and longer scales. Hence the viscous and 
higher order terms should not play a significant role in 
the observed behaviour. 

To further support the conclusion that the viscous con- 
tributions are insignificant, we will consider two ways of 
incorporating effects of these terms: (i) through the ad- 
dition of an artificial viscosity (as defined in [H, IsBj ) 
and (ii) by comparing results for higher temperatures, 
where the viscous and higher terms become less relevant 
(as described above), and showing that the observed tur- 
bulent behaviour remains unchanged both qualitatively 
and quantitatively, e.g., see figure [5l Notice that one 
could add terms up to second order in the derivative ex- 
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pansion (third order in derivatives) as suggested in j32| 
(which provides a modern extension of the method sug- 
gested in [m ) . This approach effectively amounts to con- 
trolhng and reducing gradients in the solution through a 
dynamical equation governing second order gradients of 
the flow. As we will see, the natural evolution of the 
system to longer wavelengths renders these approaches 
unnecessary. 

In what follows we describe the system of equations 
considered, discuss useful monitoring quantities — in- 
cluding the conserved quantity associated to the enstro- 
phy — as well as a brief summary of stationary solutions 
inferred from the dual gravitational picture. 



A. Evolution Equations 

The system follows from the local conservation of the 
stress-energy tensor of a perfect fluid. 



(2) 



along with the condition of conformal invariance that re- 
quires the stress tensor to be traceless, i.e., = 0. This 
latter condition fixes the equation of state as 



(3) 



for a conformal fluid in 2+1 dimensions — compare with 
([1]) for general dimensions. 

The set of dynamical variables we have chosen to evolve 
the system are: {p, i^^}, where p = log{p^^^). The 
remaining component of the three- velocity, u^, which 
can be identified with the Lorentz factor 7 {u^ = 7 = 
(1 — is obtained at each time-step through the 

normalization condition U/j^u^ = — 1. In equilibrium, the 
temperature T of the system is related with the above 
variables by p = T^, or equivalently, T = 6^.(56^ 

The system of evolution equations obtained for our dy- 
namical variables then reads, 
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dtp = —\zVp^zd-^Vz) , (4) 

dtUi = — ( [adi + UiV] p + aVui - zuid + ^Vz\b) 
a7 I 2 J 



The spatial derivatives along the flow are denoted by 
V = u^Dk^ while = DkU^ is the spatial divergence 
and Dk is the covariant derivative associated with the 



spatial metric hij. 
a = 1 + 7^. 



We also have defined, 2: = 7 and 



B. Equilibrium Configurations 

Equilibrium configurations for the toroidal case are 
simply constant temperature/constant velocity fluid 
flows and we adopt conditions as described in the ap- 
pendix. The equilibrium states for conformal fluids on 



the two-sphere can be obtained straightforwardly by 
requiring no entropy production within the first sub- 
leading (dissipative) order in the fluid expansion at equi- 
librium. This imposes a restrictive condition for the shear 
viscosity (i.e., (Tab = 0). A simple family of such config- 
urations corresponds to rigidly rotating fluids, which in 
our variables are given by 

= , ii^ = 7a;o , T = 7t (6) 

where 7 = {1—uJq sin^{9))~^^'^ is the Lorentz factor. The 
two (constant) parameters characterizing these solutions 
are: the cjq, the angular rotation rate, and r, the lo- 
cal temperature measured by co-moving observers. Note 
that implicitly, we have set the radius of the S'^ to one 
here, i.e., the proper length of the equator, = 7r/2, is 
simply 27r. 

The thermodynamics and the local stress tensor of 
these solutions has been found to be in precise agree- 
ment with thermodynamics and boundary stress tensor 
of the spinning black holes In this reference, the 

authors compare conformal fluids on spheres of arbitrary 
dimensions with large rotating black holes on AdS spaces. 
First, global thermodynamical quantities are compared. 
Then, appeahng to the duality, a comparison is made be- 
tween the local stress tensor of the fluid configuration and 
the boundary stress tensor for the most general rotating 
black hole in AdSd^rii as given in [s^- The relevant 
black hole solution for our case, the one corresponding 
to Kerr AdS 4,^ is labeled by two parameters a and r+, 
related with the angular momentum (per unit mass) and 
the horizon radius of the black hole, respectively. Perfect 
agreement was found between these two theories in the 
large r+ limit, upon the following identifications [3 61] : 



a, 
3 r+ 
47r 



(7) 
(8) 



Thus, the static configuration with cjq = is the fluid 
dual to AdS/^ Schwarzschild black hole, and the rotating 
fluid configurations (cjq 7^ 0) are dual to the AdS/^ Kerr 
geometry. 



C. Conservation Laws 

Conserved quantities are invaluable tools to analyze 
the behaviour of any given system and, as discussed, their 
existence can imply a particular behaviour of the fluid 
flow. We here discuss these quantities which we monitor 
in our efforts to understand the behaviour of the system, 
in particular the phenomenon of turbulence. In particu- 
lar, given a conserved current (i.e., V^J^ = 0), the 
existence of a conserved charge immediately follows. 



/ 



J^dH = const . 



(9) 



Hence in the following, we identify various conserved cur- 
rents and charges for the present system. 
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Given a killing vector field Kjj^ and the conserved stress 
tensor, the current = T^^Ky^ is automatically con- 
served. Thus, from the killing vectors, ^ ^ and 
C ~ we construct the first two conserved quantities 
that are identified with the total energy and angular 
momentum of the fluid, respectively: 



E = I f p(37'-l)dE, 



(10) 
(11) 



From thermodynamical considerations, one knows that 
the local entropy s = p^^^, is conserved along the flow 
direction (in a perfect fluid). Thus, the current = su^ 
is conserved and total entropy conservation follows. 



/ 



(12) 



Another quantity of interest is the vorticity which 
arises as a purely geometric property (i.e., independent 
of the equations of motion). It is constructed by tak- 
ing the exterior derivative to the flow velocity, resulting 
in the vorticity two-form u^^y = d^jj^Uy-^. In 2+1 dimen- 
sions, one can define a naturally conserved current just 
by taking its Hodge dual, i.e., we define = e'^^^uo^y. 
The total 'circulation', which remains constant through- 
out the evolution, is then 



C = I W^dT.= I . 



e^'^Uij dY. . 



(13) 



However, note that the conservation law is 'topological' 
in this case and so this conserved charge actually vanishes 
since C = — J^2 du = 0. 

Finally, we define the enstrophy, whose conservation 
plays a crucial role on two-dimensional turbulence and 
its inverse cascade scenario. This quantity, for incom- 
pressible non-relativist ic fluids, is just the integral of the 
square vorticity field. Carter [38| has generalized the con- 
cept for relativistic fluids in three spatial dimensions, and 
we extend it here for two-dimensional relativistic confor- 
mal fluids. 

Let us begin with V^T^^ = with the stress tensor 
given in ([2]) for a conformal fluid in d spacetime dimen- 
sions (i.e., v — pI {d—V). For convenience we project these 
conservation equations along and orthogonal to the flow 
velocity, 

^""d^P = -^P{^>^^'')^-^^P^. (14) 



(15) 



which is built to satisfy the Carter- Lichnerowicz equation 
of motion 



0. 



(17) 



which is equivalent to ([T5j). Note that the inclusion of 
p^/^ in (p!6|) is crucial to produce this formulation of the 
equations of motion. Then, from the Cartan Identity^ is 
straightforward to show that 



£au^^ = Au • dl^ + d(Au • 17) = . 



(18) 



Here the first term vanishes because Q, = d{p^^^u) is 
an exact form, while the second term vanishes by the 
Carter-Lichnerowicz equation (pTj) . This means that the 
two-form Q does not change along the flow direction. 

The latter observation motivates one to look for a new 
conservation law, however, we will only be able to con- 
struct the desired result by using an identity which only 
holds for d = 3, i.e., two spatial dimensions. In this case, 
we write the following current: = p~'^^^{Q'^^Qa(3)u^ • 
To establish that this current is conserved, we make use 
of the identity. 



(19) 



where ^} = Vt^yVt^^ . The latter holds for two spatial 
dimensions since we can write, Q.^y = e^y^l^, for some 
arbitrary vector /'^. The condition Cl^yU^ = then re- 
quires that be proportional to u^, i.e., = (_u^. One 
can then find the proportionality factor ^ in terms of 
to show (plQl) holds. With this relation, we have 



= 2p- 
= 2p- 



i {en^ -2n''°'n^p{Vau^)} =o. 



Here we have used (fHj) in the second line and (fTSj) in 
the third hne, while the final vanishing follows from ()19p . 
Therefore, the enstrophy, 



> Z 



(20) 

is conserved, where the first term in the expression is 
just the vorticity two- form squared (as expected from 
the non-relativistic version), and the second one, which 
involves the acceleration = u^VyU^^ accounts for 
the fact that the fluid world-lines are not necessarily 
geodesies. One important point to emphasize from (|2Q|) 
is that the factors of p cancel out in the end, so that the 
final expression for the enstrophy is independent of the 
energy density. 



where O is the (full covariant) divergence; a^, the acceler- 
ation; and P^^ = g^^ -\-u^u^ ^ the projector perpendicular 
to u^. Next, consider the two- form. 



(16) 



D. Numerical approach and setup 

Our goal is to study the dynamical behaviour of a per- 
turbed, otherwise stationary, fluid configuration dual to 
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the Kerr/Schwarzchild black hole in global AdS^ or the 
Poincare patch. We next describe the different compo- 
nents of our numerical implementation. 



1. Initial Data 

Initial data for perturbations of the stationary fluid 
configurations are directly added to the fluid velocity 
such that eddys are induced. This is straightforwardly 
achieved by considering perturbations which, in the co- 
moving frame of the fluid, describe space varying veloc- 
ities that change sign in a smooth manner along some 
particular direction chosen. In the case of the torus, this 
is straightforward as described in the appendix. For the 
spherical case, we adopt a perturbation in the fluid ve- 
locity as: 



(21) 



for some small value 6 and a function cUp defined on the 
sphere (generally chosen to be one of the spherical har- 
monic basis functions Y^{0^(j))). In this expression, 7 
is the local Lorentz factor defined for the full angular 
velocity cj, including the Up contribution. For future ref- 
erence, it is useful to indicate the initial vorticity density 
associated to this flow. 



{multi-block approach). To ensure a correct transfer of 
information among the different grids we follow the tech- 
nique described in [39], which relies on the addition of 
suitable penalty terms to the evolution equations to pre- 
serve the energy norm through the whole sphere. This 
technique 'penalize' possible mismatches between values 
the characteristic fields take at interfaces and enforce con- 
sistency through suitably introduced driving terms. Here 
we follow the strategy introduced in 40 1 
tension of the method introduced in 41, 



which is an ex- 

Si. 




FIG. 1: Cubed Sphere Coordinates. A total of six Cartesian 
patches are employed to cover the sphere. Only patch bound- 
aries coincide at common points. 



[sm\0)j{uJo^Sujp{0,^))] 



Vh sm(6>) 
= (cjo + (^cjp)7^cos(6>)(2 -cj^sin^((9)) 
+(^7^ sm{e)deujp . 



From this expression, we see that the rigid rotation com- 
ponent (|6j) of the flow introduces a Yi component to 
the vorticity at order ujq and higher £ contributions also 
appear at higher orders in ujq. Thus for sufficiently fast 
flows, the background contribution renders analyzing tur- 
bulent behaviour through vorticity more delicate. 



Introducing coordinates {t^x^y} with (x^y) to label 
points in each Cartesian patch, the metric in each one 
reads. 



dt^^^ {(1 + y^)dx^ + (1 + x^)dy^ - 2xydxdy] 



where D = 1 ^ x^ ^ y^ . The non- vanishing Christoffel's 
symbols are given by. 



-2^ 



]J ^ xy jj yy ]J ^ xy jj 



2. Grid Scheme 

We now discuss details of our grid implementation and 
results. Since the qualitative behaviour is the same in 
both topologies considered -but the implementation is 
more involved in the spherical case- from now on we 
concentrate on the spherical flows to simplify the pre- 
sentation and defer details of the toroidal case to the 
appendix. 

Since the topology of our computational domain is 6*^, 
we employ multiple patches to cover it in a smooth way. 
A convenient set of patches is defined by the cubed sphere 
coordinates. There are six patches with coordinates pro- 
jected from the sphere, and each of this patches consti- 
tute a uniform grid — see figure [1]). These grids are 
defined in a way such that there is no overlap and only 
grid points at boundaries are common to different grids 



3. Numerical Scheme and Stability 

In order to construct stable finite difference schemes for 
our initial value problem we use the method of lines Q . 
This means that we first discretize the spatial derivatives 
(constructing some suitable finite difference operators) 
and obtain a system of ordinary differential equations 
for the grid functions. To ensure the stability of the nu- 
merical scheme we use the energy method described on 
(35} . We employ (fourth/second-order accurate at inte- 
rior/boundary points) finite difference operators satisfy- 
ing summation by parts (the discrete analogue of integra- 
tion by parts) and deal with interface boundaries with 
appropriate penalty terms at the interfaces as described 
in [39|. For the time integration, we adopt a 4^^ order 
Runge-Kutta algorithm. 

All simulations were performed using 81x81 grids for 
each of the six patches, giving a total number of around 
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40.000 grid points to cover the entire sphere. We have 
confirmed convergence through several tests, in partic- 
ular we adopted different initial data, and studied the 
numerical solutions with increased resolution with each 
grid having 321x321, 161x161, 81x81 points (labeled by 
4,2,1 respectively). With each obtained solution, denoted 
by U^\t,e,(t)) (for i = 4,2,1), we calculated the conver- 
gence rate p as, Q{t) = ||^(2)_^(i)||^^ ^ 2^- For the cases 
considered, the obtained rate was in very good agreement 
with the expected for resulting 3^^ order implementation 
constructed. 



III. RESULTS 

We analyze the dynamics of small perturbations 
around the stationary fluid configurations duals to 
Schwarzschild and Kerr geometries on AdS^. We sep- 
arately study these two cases, starting with a qualita- 
tive description of the system evolutions, for generic per- 
turbations. We show in figures 2] and [6] the sequence 
on the evolution of the vorticity field for the two cases. 
Then, in order to gain some insight into the turbulent 
behaviour and to capture the possible cascading phenom- 
ena, we perform a spherical harmonic decomposition (up 
to i=12) of the relevant fields and compute their associ- 
ated power spectrum as a function of time. 

This signal processing analysis is derived simply from 
a generalization of ParsevaVs theorem^ which states that 
the total power of a function / defined on the unit sphere 
is related to its spectral coefficients by, 

. oo 

- / |/(^,</>)pdE = ^C^(£), 



£=0 



where are the coefficients of the expansion of / in 
the spherical harmonics. In figures 2] and [51 we plot the 
coefficients C{i) as a function of time, to analyze how 
the different modes behave during the evolution. 



A. Non-rotating case: perturbations to 
Schwarzschild 



One can basically classify the dynamical behaviour of 
the system into four different stages. In figure [21 we plot 
the vorticity field of a representative example for each one 
of them. A first stage corresponds to an initial transient 
period when the initial configuration seemingly remains 
unchanged for some time interval. (This interval depends 
on the perturbation considered, being shorter for larger 
-^'s perturbations at a fixed value of (5, but it does not 



depend on the system's temperature). Closer inspec- 
tion however reveals that non-trivial dynamics begins 
to manifest and an exponential growth of some modes 
sets in (see fig. [1]). Notice that the initial growth rate 
of these modes is approximately the same though the 
higher ^'s grow to larger values earlier. Since truncation 
errors feed higher frequencies this behaviour is not sur- 
prising. As these modes become large enough (roughly 
commensurate to the magnitude of the initial perturba- 
tion), the original symmetry of the system gets broken 
and a number of eddies arise and move around. Such an 
inst ant m ight be regarded as the beginning of turbulence 
(fig. |2(b)[ ). As the dynamics continues, the eddies gradu- 
ally turn into individual vortices and exhibit a seemingly 
chaotic motion, during a stage that we will refer as fully 
developed turbulence (fig. 2(c)[ ). With the vortices prop- 
agating around the sphere, encounters of same-sign vor- 
tices lead to increasingly larger vortex structures. During 
this stage, governed by non-linear effects, energy trans- 
fers from the higher / modes to lower ones. The process 
continues until four vortices, two of each sign, are formed 
(fig- Hd)]). 

The end-state is found to be qualitatively the same in 
all the cases we have considered, regardless of the ini- 
tial perturbation introduced and of the system's temper- 
ature. We have explored generic perturbations in the 
velocity initial data, as described on (|2T]) : from different 
single- mode perturbations (i.e., Up ~ Yf^{6^(j)) for par- 
ticular i and m), to a random combination of all of them 
up to £ = 12. We have further considered the inclusion 
of a random forcing term in the evolution equations, by 
adding a term to the right hand side of the stress 
conservation equation as: V^T^^ = This force is 
given by a random field, both in time and space (and 
thus, acting on very short scales). Even in this forced 
case, the dynamical character and properties of the so- 
lution remain qualitatively unchanged. This (in addition 
to convergence studies) indicates that our grid structure 
is playing no significant role on the obtained behaviour 
and on the final configuration attained. 

Main features of the long-term behaviour of the solu- 
tion are illustrated in figure [3l In particular, the men- 
tioned four dominant vortices of the resulting configura- 
tion. It can be noted that the temperature and energy 
densities attain local minimums at the vortices' locations 
and that the enstrophy is almost exclusively contained 
within them. 

In figure [H we have displayed a few representative 
modes of the vorticity power spectrum for an initial per- 
turbation ujp{e, (j)) = Y^Q{e, (I)), S = 0.2 and r = 100. The 
coefficients C{i) of the spectrum are normalized with 
respect to the total power and plotted (in logarithmic 
scale) as a function of time. In the figure, the modes 
£ = 9 and £ = 11 dominate the spectrum at t = — 
as expected since the vorticity involves a derivative of 
the three-velocity and the initial data for the latter has 
a single £ = 10 mode. The rest of the modes exhibit 
an exponential growth on the early stages of the evolu- 




(d) t = 950 



FIG. 2: Evolution of the vorticity field for a perturbation 
ijOp{0^(j)) — ¥10(^,0) and 5 — 0.2, on Schwarzschild (cjq). (a) 
Initial config., (b) beginning of turbulence, (c) fully developed 
turbulent stage, (d) final state. 




(a) Vorticity 
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(d)Enstrophy 



FIG. 3: Late times configuration for distinct relevant fields 
in the non-rotating case {ujq = and T ^ 100, at t = 1200). 
Notice that vortices correspond to minima in the energy. This 
behaviour is expected as stable vortex structure require a 
larger surrounding pressure that prevents its dispersion. 



time 



FIG. 4: Relevant modes in the power spectrum of the vorticity 
field for a ujp{0^(j)) = Yio{0,(t)) and S = 0.2 perturbation, in 
the non-rotating case (at temperature T ^ 100). 



the observed turbulent phenomena, and the conclusions 
we can draw from this study, might be significantly af- 
fected when such dissipative terms are taken into ac- 
count. However, recall that in the discussion around ([T]), 
we argued the perfect fluid equations will give a good 
description when LT 1. We verified the accuracy of 
this statement by comparing the system's dynamics at a 
variety of temperatures — recall that the radius of the 
sphere is fixed to be one — while keeping the initial data 

for fixed. Figure [5] illustrates the behaviour of the L2 

1/2 

norm of ^(^ i.e., [Jg^ \u^\'^dT] ' . The latter is a useful 
proxy for the solution's turbulent behaviour since in the 
absence of turbulence, it would remain zero by symmetry 
considerations. As is evident in this figure, the system's 
behaviour does not change as the temperature increases. 
In particular, the growth rate and both onset and full 
development of turbulence are essentially the same in all 
cases considered. 



B. Rotating case: perturbations to Kerr 



tion, as previously described. As the turbulent cascade 
begin, and after a short stage in which the mode struc- 
ture is rather complex, the original high-^ modes decrease 
while the lower ones increase and gradually dominate the 
flow. Particularly the £ = 2 mode in the spectrum, which 
represents the quadrupole contribution, is the dominant 
mode at late times. 

Note that the vortices in figure [3l i.e., the £ = 2 mode 
above, are essentially quasi- stationary as we have ne- 
glected here the viscous (and higher) contributions which 
would certainly affect the very long-time behaviour of the 
solution. As we noted before, one might wonder whether 



The qualitative behaviour of the system on the rotat- 
ing scenario is illustrated in figure [6l and it happens to be 
very similar to the presented above for the Schwarzschild 
case. However, we emphasize that these rotating solu- 
tions are not simply static solutions, i.e., ujq = 0, in a 
rotating frame. Recall that the rigid rotation introduces 
a background vorticity field, which one finds completely 
dominates the long-term behaviour. A key difference is 
evident immediately after the vortices are formed: the 
background rotation drags them into the main rotating 
stream and gradually separates them according to the di- 
rection of their rotation. If the background rotation flows 
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FIG. 5: Dependence with temperature: logarithm of the 
L2 norm of for different temperatures, in the non-rotating 
case. 



left to right as in fig. 6(c)[ clockwise rotating vortices ac- 
cumulate in the southern (lower) hemisphere while the 
vortices with counter-clockwise rotation migrate towards 
the opposite pole. Then, the system undergoes a merg- 
ing process (of co-rotating vortices) as in the non-rotating 
case, but now just one vortex of each sign remains at late 
times, as in fig. 6(d)[ The initial northwards/southwards 
propagation of counterclockwise/clockwise vortices can 
be also explained in terms of the equal-sign mergers of 
vortices by regarding each of the smaller vortices gen- 
erated in the flow as interacting with two much larger 
vortices induced at the north and south pole by the rigid 
rotation component of the flow. Note that the 'finaV 
two vortices are not necessarily in precise alignment with 
the rotation axis, being often the case that they remain 
orbiting around the poles. 

Main features of the long-term behaviour of the so- 
lutions in the rotating scenario are illustrated in figure 
[71 In particular the mentioned two dominant vortices 
as they oscillate around the poles, with the temperature 
and energy attaining local minima at the vortices' loca- 
tions and being concentrated near the equator where the 
fluid velocity is larger. The enstrophy is again almost 
exclusively contained within the vortices. 

In figure [8] we display representative modes of the vor- 
ticity spectra for two different parameters of the rotat- 
ing solution with r = 100 and a perturbation set by 
ujp{e,(l)) = y?o(<9,0) and S = 0.08. Here, the ^ = 1 
mode represents the background rotation contribution 
while the i = 9 and i = 11 modes are the ones associ- 
ated with the perturbation, as discussed previously. We 
should recall here that the coefficients in the spectrum 
are normalized to unity, and thus it becomes clear from 
the plots that the long term dynamics is being completely 
dominated by the rotation. 

For the case ujq = 0.1, the configuration starts with 



FIG. 6: Evolution of the vorticity field for a perturbation 
uOpiO^cj)) — YiQiO^cj)) and 5 — 0.08, on a rigid rotation with 
cjo = 0.1. (a) Initial config., (b) beginning of turbulence, (c) 
fully developed turbulent stage, (d) final state. 




(a) Vorticity 



(b) Temperature 




(c) Energy 
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FIG. 7: The late time configuration for distinct relevant fields 
in the rotating case (cjo = 0.1 and T ^ 100, at t = 3600). 



a perturbation comparable in magnitude with the back- 
ground motion. As the vortices form and turbulence be- 
gins, the high-^'s modes decrease displaying a cascading 
phenomena into lower I modes. In particular this cas- 
cade progresses towards the ^ = 1 mode that eventually 
governs the state of the system. On top of this rotation 
dominated flow, the two vortical structures mentioned 
above persist and are represented in the spectrum by a 
combination of higher modes, predominantly the 1 = 2 
and ^ = 4 modes. 

Interestingly, for faster initial background rotations 
(i.e., ujQ > 0.5) while keeping the perturbation ampli- 
tude S fixed, the flow is already dominated by the i = 1 
mode from the beginning and the turbulent stage takes 
longer to appear and fully develop. It seems that a strong 
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FIG. 8: Relevant modes in the power spectrum of the vorticity 
field in the rotating case (at temperature T ^ 100), for two 
rotation parameters ojq. The initial perturbation is given by 
ujp{0, 0) = Y^o{0, 0) and 6 = 0.08. 



rotating stream in the background hinders to a certain 
degree the formation and merging of vortices. 



IV. DISCUSSION AND FINAL WORDS 

In this work we have analyzed relativistic conformal 
fluid flows on both S'^ and backgrounds and found 
first that turbulence naturally arise in these flows and 
second that it gives rise to a 'inverse' cascade from 
shorter to longer wavelengths. These results not only ex- 
tend common observations in Newtonian hydrodynam- 
ics but also have tantalizing implications for the be- 
haviour of gravity in four-dimensional AdS spacetimes. 
Indeed, the AdS/CFT correspondence implies that for 
sufficiently high temperatures/length scales, the confor- 
mal fluid flows studied here have a dual description in 
terms of gravitational perturbations on Schwarzschild 
or Kerr black holes in AdS4, either with a spherical 
or planar horizon. The turbulent behaviour observed 



here implies, through the duality, that gravitational per- 
turbations in this limit should cascade to smaller fre- 
quencies. This constitutes a prediction obtained within 
holography which had not been previously anticipated 
on flrm grounds [57j. Indeed, previously identifled in- 
stabilities in AdS, super- radiance (e.g., [43|) and 'weak- 
turbulence^ imply a frequency shift towards higher 
frequencies. Therefore the cascading phenomena ob- 
served here implies an altogether new behaviour on the 
gravitational side. Furthermore, as described in [sl, [H, [36[ , 
the full metric of the corresponding spacetime can be ob- 
tained and the implications of this cascading behaviour 
analyzed through suitable geometric quantities. We de- 
fer such tasks to a forthcoming work. We close noting 
that this possible cascade behaviour, flrst raised in [5[ 
and demonstrated here (and presented in e.g.. 1461 ) has 
now been observed in the Poincare patch case [43] . 

Of course, another fleld theory 'prediction' is that large 
AdS black holes in five and higher dimensions will also ex- 
hibit turbulence but this chaotic behaviour will give rise 
to a 'standard' cascade to shorter wavelengths in these 
cases. From a gravitational perspective, this apparently 
generic behaviour for 'hot' horizons in AdS is completely 
unexpected and calls for a better understanding within 
gravity itself. 

A step towards understanding the distinction between 
gravity in four and higher dimensions can be obtained by 
recalling that enstrophy conservation is what drives the 
hydrodynamics in 2+1 dimensions to exhibit the inverse 
cascade. One can then exploit the duality to translate the 
enstrophy into geometrical variables and to understand 
the implications of its conservation on the gravitational 
side. That is, the conservation of enstrophy in the fiuid 
description implies a quasi-conserved quantity exists in 
the bulk gravity theory. Further, as we have seen here, 
the system displays a rich dynamical vortex configura- 
tion that merge towards a long lived state described by 
relatively few, long wavelength vortices. Isolated vortices 
in conformal fiuids in 2+1 dimensions have been studied 
in and their properties identified. Again using the 
fiuid/gravity correspondence, one could produce a grav- 
itational description of these quasi-stationary vortices. 
That the understanding of these or other geometrical 
quantities might help shed new light in turbulence phe- 
nomena is definitively an exciting prospect (see e.g., |49|). 

Certainly, the full gravitational description will nat- 
urally incorporate dissipative contributions in the fiuid 
fiows. While we argued and quantitatively demonstrated 
that this dissipation will not modify the essential features 
of the turbulent behaviour at sufficiently high tempera- 
tures. In Newtonian hydrodynamics, the onset of turbu- 
lence is discussed in terms of the Reynolds number Re. 
The typical benchmark for the onset of turbulent fiows 
is that Re have value of a few thousand. In relativistic 
hydrodynamics, the latter may be estimated as [50| 



TL 

Rer = 47rTL, 

r]/s 



(22) 
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where L is a characteristic length scale in the flow. In the 
last expression above, we have substituted the celebrated 
holographic value r^/s = l/47r [11,0, Uli — we would have 
rj/s — j3 for the general conformal fluid in ([T]). As this 
expression illustrates, we can produce arbitrarily large 
values of Re by increasing T (while keeping the system 
size flxed). Hence to observe turbulent behaviour, we are 
again naturally pushed to the regime where the hydro- 
dynamic gradient expansion works well and our perfect 
fluid model becomes a good approximation. Of course, 
the viscous effects will deflnitely modify the very long- 
time behaviour observed here. For example, the conser- 
vation of the enstrophy is only true to flrst order in the 
hydrodynamic gradient expansion. 

Of course, another interesting extension of the present 
investigation would be studying further the details of the 
turbulent cascade, in particular, the Kolmogorov scal- 
ing exponents. Figure [13] illustrates some preliminary re- 
sults, which suggest that the 'Newtonian' kinetic energy 
will scale with the expected exponent of —5/3. However, 
various caveats must be noted. First, a clean easy-to- 
distinguish Kolmogorov-type reasoning applies to a sit- 
uation where the turbulent flow is driven by an exter- 
nal force at some high frequency (in the present case 
of 2+1 dimensions) and viscous dissipation also damps 
the energy flow on much longer time scales. This sce- 
nario must be contrasted with the freely-decaying tur- 
bulence (i.e., without any driving force) which flgure [T3l 
describes. In freely-decaying turbulence determining the 
intertial range is more delicate but in such range (which 
shrinks as time proceeds) the —5/3 slope can be distin- 
guished. Additionally, Kolmogorov's scaling arguments 
are made in a Newtonian context and can at best be re- 
garded as approximate for the relativistic fluids studied 
here. Fortunately, the appropriate exact scaling relations 
applicable for relativistic hydrodynamic turbulence have 
been derived in [50]. The present simulations provide 
a framework for further studying these relativistic rela- 
tions. More generally, however, the most exciting possi- 
bility would again be if a holographic perspective could 
provide new insights into the issues surrounding these 
turbulent cascades. 
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Appendix A. 

For comparison purposes, we also consider scenarios 
related to the dual case of an AdS^ black brane solu- 
tion where the non-radial directions are compactifled on 
a torus of size D. In Feffer man- Graham coordinates, la- 
beled by (t, X, ?/), the asymptotic metric is simply the flat 
one and the equations of motion can be straightforwardly 
implemented. 



A. Stationary solutions and initial data 

Equilibrium conflgurations are given simply by con- 
stant flows at a given (arbitrary) temperature T. For 
simplicity we consider a torus with domain [0,D]^. We 
then restrict to initial conflgurations with a flow along 
X, i.e., = 5^uq^ and introduce generic perturbations 
8u^ to this flow. For concreteness, we describe here two 
particular cases where the initial three-velocity is given 
by: case (A) with a perturbation of compact support, 

= 51 {u^ + 5u sin(2^^/3) {y - Lf) (A-1) 

for < ^ < L and = S^uq otherwise; and case (B) 

= {uo + Su sm{167ry/D)) . (A-2) 

We note however the qualitative features observed in all 
cases considered remain the same. 



B. Turbulence, cascading behaviour and 
temperature dependence 

We performed a series of numerical experiments adopt- 
ing Uo = 0, 0.1 and 0.5 and 6u = 0.01, ••• , 0.05 (i.e., 
perturbations from 2% to 50%). and considered temper- 
atures T G [1, 10^] with L = 10 and typical grid sizes of 
[401,401] (though consistency of behaviour was checked 
with grids 1.5 and 2 times better resolved). 

As in the case of the sphere, the dynamics display a 
turbulent behaviour leading to the development of large 
vortices. As an illustrative example. Figure [9] displays the 
vorticity of the system at different times for case (A) with 
Uo = 0.1^ du = 0.01 and T = 1. Early on, the behaviour 
is seemingly stationary but as turbulence develops the 
initial symmetry is completely broken and vortices arise 
which grow as they merger leading to a configuration 
described by long wavelengths. 

As we have done in the S'^ case we also study the sys- 
tem's behaviour upon variation of temperature to ensure 
terms neglected in our study do not significantly affect 
the dynamics obtained. As in the previous case, a use- 
ful proxy to monitor the evolution is the (L2) norm of 
u^ which, in the absence of turbulent dynamics, would 
remain zero. Figure [10] illustrates that the observed be- 
haviour is essentially unchanged with temperature. 
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FIG. 9: Representative snapshots of the vorticity. As time 
progresses the initial configuration is strongly disturbed by 
the formation of vortices. As the dynamics continues, larger 
vortices are formed. 
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FIG. 10: Logarithm of the L2 norm of for different tem- 
peratures. Essentially the behaviour observed is unchanged 
with temperature. 



Last, we also considered varying the size of the torus 
by increasing D by factors of 1.5 and 2 in case (A) and 



observed the same cascading behaviour. 

To study how the energy cascades as time progresses, 
we compute the Fourier transform of energy density 
T^^ and the vorticity density . Since our computa- 
tional grid is discretized by points {x\^x^j} (with i^j = 
I,-- - ,A^), we denote n = {i^j}^ and write T^^(x) = 
T^^{n) and the vorticity W^{x) = W^{n). The Fourier 
transform of these quantities is given by 

f°°{K) = ^e-2-^-^T™(n), (A-3) 

n 

W°{K) = ^e-2"-^-^T^°(n). (A-4) 

n 

Here K = {i^i,i^2} and i^i,i^2 e [0,7V - 1]. Clearly, 
T^^{K) and W^{K) are also represented on a TV x TV grid. 
It is convenient to re-express the transformed quantities 
as functions of K = \K\. From the functions T^^{K) and 
W^{K), we define T^^{K) and W^{K) as follows 

T'\K) = Y^\T'\K,)\, (A-5) 

i 

W°iK) = J2\W°iKi)\, (A-6) 

where the sums run over Ki^s satisfying K < \Ki\ < 
K With this approach an effective wave-number grid 
of length ^/2N is obtained. Figure [TT] illustrates a clear 
cascade to lower wave- numbers as time proceeds. In the 
case of T^^{K) the dominant mode is given by = 1 
while for the vorticity both K = 1 and K = 2 are of 
similar magnitude. As in the S'^ case, this is expected as 
the vorticity field is obtained by taking a single derivative 
which increases the mode content by one. 

It is interesting to monitor whether the observed be- 
haviour is consistent with the 'standard' Kolmogorov ex- 
pectation. To illustrate this, we compute the Fourier 
transform of the 'Newtonian' kinetic energy per unit mass 
E l/2v^. To make sense of this limit within our rel- 
ativistic description, we chose uq = and Su = 0.03, 
whose solution is such that \v\ is bounded by 0.04 for 
all times, thus the fluid's motion stays far from relativis- 
tic speeds. Figure [12] shows the vorticity behaviour for 
this case and [T3l illustrates the Fourier transformation of 
this kinetic energy for representative times. As time pro- 
ceeds the energy in higher frequency modes diminishes 
while the opposite behaviour is observed for the lower 
ones. Additionally, at inter midiate frequencies, the en- 
ergy exhibits a behaviour with frequency consistent with 
a slope of —5/3, the expected Kolmogorov exponent. 
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